A Model for Ferromagnetic Nanograins with Discrete Electronic States 
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We propose a simple phenomenological model for an ultrasmall ferromagnetic grain, formulated 
in terms of the grain's discrete energy levels. We compare the model's predictions with recent 
measurements of the discrete tunneling spectrum through such a grain. The model can qualitatively 
account for the observed features if we assume (i) that the anisotropy energy varies among different 
eigenstates of one grain, and (ii) that nonequilibrium spin accumulation occurs. 



PACS numbers: 73.23.Hk, 75.50.Cc, 73.40. Gk 

What are the properties of individual quantum states 
in the electronic excitation spectrum of a nanometer- 
scale ferromagnetic particle? This is becoming an in- 
creasingly important question, since the size of mem- 
ory elements in magnetic storage technologies is decreas- 
ing extremely rapidly and particles as small as 4 
nm are coming under investigation ||. In this size 
regime, the excitation spectrum becomes discrete; in- 
deed, Gueron, Deshmukh, Myers and Ralph (GDMR) 
H have recently succeeded to resolve individual quan- 
tum states in the spectrum of ferromagnetic Cobalt 
nanograins, using single-electron tunneling spectroscopy. 
They found complex nonmonotonic and hysteretic energy 
level shifts in an applied magnetic field and an unex- 
pected abundance of low-energy excitations, which could 
not be fully understood within the simple models used 
for ferromagnetic nanograins so far [0,0. 

In this Letter, we propose a phenomenological model 
for ferromagnetic nanograins that is explicitly formulated 
in terms of the discrete states occupied by the itiner- 
ant conduction electrons and capable of qualitatively ex- 
plaining the observed features. The model is similar 
in spirit to that advanced independently by Canali and 
MacDonald but our analysis includes two further in- 
gredients beyond theirs: (i) mesoscopic fluctuations of 
the anisotropy energy (i.e. it may vary among different 
eigenstates), and (ii) nonequilibrium spin accumulation. 

Experimental Results. — GDMR studied Co-particles 
1-4 nm in diameter. Assuming a hemispherical shape, 
the number of atoms in such grains is in the range 
N a w 20-1500, and the total spin, s « 0.83AT a §], thus is 
so ~ 17— 1250. In GDMR's devices, a grain is connected 
to two aluminum electrodes via aluminum oxide barriers. 
Its tunneling conductance consists of a series of distinct 
peaks (see Fig. 2 in ||), whose positions yield a set of tun- 
neling energies of the form Q AE^ = E^ 1 — Ef each 
corresponding to the energy cost of some rate-limiting 
electron tunneling process \i) N — > \f) N±1 onto or off 
the grain. Here \i) N denotes a discrete eigenstate, with 
eigenenergy , of a grain with N electrons, etc. 

As the magnetic field is swept, the resonances for Co- 



grains undergo energy shifts and crossings (Fig. 3 in ||, 
0). The resulting tunneling spectra have several proper- 
ties that differ strikingly from those of previously-studied 
nonmagnetic Al and Au grains |8|||: 
(PI): Many more low- energy excitations are observed 
than expected: For all values of magnetic field, the mean 
level spacing is c? b s ^ 0.2 meV. This is much smaller 
than expected from the naive single-particle estimates 
4.6 eV/so or J min « 1.2 eV/s (with s < 1250) 



for the majority- and minority-spin mean level spacings 
near the Fermi energy of Co Jl(]|| . 

(P2): In the small-field regime ((j,qH < 0.2 T), discon- 
tinuous hysteretic switching occurs at a certain switching 
field H sw , due to a sudden change in direction (henceforth 
called "reversal") of the magnetic moment. Moreover, 
the iJ-dependence of tunneling resonance energies has 
continuous non-monotonic variations, which differ seem- 
ingly randomly from level to level (Fig. 3 in ||, 0). 
(P3): In the large-field regime (|-Ff| 3> \H SW \), the res- 
onance energies depend roughly linearly on H, with H- 
slopes that almost all have the same sign for a given 
grain; in particular, slopes of opposite signs due to 
Zeeman-splitting of spin-up and spin-down levels ||[| are 
not observed (Fig. 4 in ||, [§). 

Point (P2) indicates immediately that an independent- 
electron approach to the energy levels is not sufficient, 
because the energy of a given state depends on the ori- 
entation of the magnetic moment produced by all the 
electrons within the particle. We shall argue that points 
(PI) and (P3), too, are related to the many-electron spin 
structure within the particle. 

Model Hamiltonian. — We propose to model a 
nanoscale magnet with discrete excitations by the fol- 
lowing "minimal" Hamiltonian: 7i = H.q + Tic + Tt cxc h + 
^Zcc + Wanisj where Tic is the Coulomb charging energy 
for a nanoparticle containing N electrons, and 



u 



s-s 



(1) 



^Zoo = —hS z 



££sf/c^, ( 2 ) 

ab ij 



l 



with h = g c fffiBlM)H . Here Ho describes the kinetic en- 
ergy of a single band of single-electron states \j, a) , la- 
beled by a discrete index j and a spin index a — (1,1), 
with the spin quantization axis chosen in the z-direction. 
The exchange, Zeeman, and anisotropy terms, 7i CX ch, 
Hzcc and H an i S , are functions of the level- j spin oper- 
ators Sj = 7;J2a'a c 'ja (T cra' c jv' ( wnere °~ a are Pauli ma- 
trices, with a = x,y,z), so that S — ■ Sj is the total 
spin vector. Tt cxc h is a rotationally invariant term which 
models the effects of an exchange field and forces the 
system to adopt a non-zero total ground state spin, say 
So- On account of this term, spins aligned parallel or an- 
tiparallel to (S) may be thought of as forming "majority" 
and "minority" bands, which effectively rotate rigidly to- 
gether with the magnetization direction. We shall take 
the mean level spacings (d ma j, <i m m) near the respective 
Fermi energies, and the exchange-splitting of the Fermi 
energies, Af = £F,maj — £F,min (~ 2 eV for Co), as char- 
acteristic parameters of the model. The magnitude of 
U may then be estimated as U w Ap/sg, since stabil- 
ity of the ground state spin sq implies B the relation 
Ap = U(s + 1/2) + d , where d (~ 1/so) is a small, 
grain-dependent energy satisfying — (d ma j —U/2) < d Q < 
dm™ — U/2. Hzcc describes the spin Zeeman energy in 
an external magnetic field H = Hz. Finally, 7Y a nis mod- 
els the combined effects of crystalline, shape and sur- 
face anisotropies, etc, in terms of a hermitian, traceless 
tensor Ko/CK" = 0], which describes the energy cost 
for rotating the various spins Sj. We split the tensor 
into an "average" and a "fluctuating" part by writing 
K ab _ R ab + k ab_ The Ra b part dominates since all lev- 
els contribute coherently, and, assuming K ab oc 1/Vol, it 
makes an extensive (oc Vol) contribution to the total en- 
ergy. The simplest non-trivial form that this term might 
take is a uniaxial anisotropy 

Wuni = ~k N (S -n^/so, (3) 

where h is the unit vector in the easy-axis direction 
(at, say, an angle 9 from z) and fcjv(> 0) is a volume- 
independent constant. The fluctuating term kf- causes 
the total anisotropy energy to depend on which single- 
particle levels are occupied. It is a new ingredient rel- 
ative to previous models of magnetic switching, which 
required only a single anisotropy energy function for the 
whole system, as is appropriate when only the ground 
state magnetic properties are pertinent . 

Basis states. — It is convenient to use the eigenstates 
of H(fC = 0) to construct a set of "bare" basis states. 
Since [H, S] = 0, these states can be grouped into spin 
multiplcts that are labeled by their S ■ S and S z eigenval- 
ues, say s(s + 1) and m. For example, the bare ground 
state of H(K, = 0) for given N, s and h (> 0), say 

\a,8)Z = ncj T n c !ii vac >' ( 4 ) 

3=1 3=1 



is a member of a spin multiplet of 2s + 1 states, \s, m) oc 
(5_)(«-fn) \ S) s)o, Here n yi = N/2 ± s, and S- = 
S x — iS y is the spin-lowering operator. For K ab ^ (but 
still kfj = 0) the true low-energy eigenstates, \s,m) , 
are linear superpositions of the bare states in the multi- 
plet s, m)g (with |s, m) N — > |s, m)^ as \K ab \/h — > + ). 
We shall call the states |s, m) the spin wave multiplet, 
since each can be viewed as a homogeneous spin wave. 
By creating additional single-particle excitations, other, 
higher-energy multiplets can be built. However, their 
eigenenergies lie higher than those of the spin wave mul- 
tiplet \s, m) N by an amount which is at least of order the 
single-electron level spacing, i.e. rather large compared 
to d bs [cf. (PI)]; thus the mechanism causing the ob- 
served abundance of low-energy excitations, whatever it 
is, must have its origin in spin excitations, not purely in 
single-particle excitations. 

Anisotropy fluctuations. — Let us first turn to the be- 
havior of the tunneling resonances for small magnetic 
fields [see (P2)]. The jumps at the switching field have 
been attributed to a sudden reversal of the nanoparticle's 
magnetic moment [p| , which occurs when the energy bar- 
rier between a metastable state and the true ground state 
is tuned to zero by the applied field. To illustrate how 
such jumps arise in our model, let us (for simplicity) take 
the anisotropy to be uniaxial [Eq. (gj)] and consider the 
case in which the changing magnetic field only rotates the 
total spin moment, without changing its magnitude [ jl2| . 
We have numerically diagonalized 7izco + W un i as func- 
tion ofh/kN for Si = 1000 and st = Sj±l/2 to determine 
the groundstate energies and the corresponding tunneling 
energies AE^ (Fig. [I]). The latter indeed show a jump 
at h S w However, if we neglect anisotropy fluctuations 
by choosing fcjy = fcjy+i (Fig. |, solid lines), the AEf. 
lines also have two unsatisfactory features: (i) An upward 
(downward) jump in AE^h) as \h\ increases past \h sw \ 
is always followed by a positive (negative) large-ft. slope, 
whereas it is observed experimentally (e.g. Fig. 3(a) in 
H) that either upward or downward jumps can occur 
for states having a given large-/i slope; and (ii), beyond 
the switching field, the dependence on h is monotonic 
(close to linear), in disagreement with recent data (P3) 
M. All attempts we made to explain such behavior by 
choices of K ab corresponding to more complicated than 
uniaxial anisotropies, or by higher order terms such as 
K abcd S a S b S c S d || , were unsuccessful. 

Now, the very fact that the field-dependence of each 
resonance in HQ differs so strikingly from that of all 
others implies that the anisotropy energy fluctuates sig- 
nificantly from eigenstate to eigenstate, which we asso- 
ciate with kf b 7^ in our model. Although the con- 
tribution of such random fluctuations to the total en- 
ergy is non-extensive, their effect on energy differences in 
which extensive contributions largely cancel, can be very 
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significant. A detailed statistical analysis of anisotropy 
fluctuations is beyond the scope of this paper. Instead, 
we shall mimic the effects of kf^ ^ by simply using 
two different anisotropy constants in W un i, say fc/v and 
fcjv±i = fcjv + Sk±, for N- or (N ± l)-electron states, fcjv 
can be estimated from the switching field fc^ w ^q^bH sw 
[cf. Fig.[l[] yielding fc^ « O.OlmcV @. Now, as illus- 
trated in Fig. [l], 5k±/k^ in the range of a few percent 
is sufficient to reverse the sign of the energy jumps at 
Hsw- Note that Sk± ^ also causes the spectral lines 
to exhibit rather strong non-monotonic "kinks" near h sw , 
whose amplitudes are of order So5k±. Qualitatively sim- 
ilar non-monotonicities have indeed been observed re- 
cently 0, with kink amplitudes on the scale of a few 
0.1 meV, in rough agreement with so5k± for sq ~ 1000. 

Anisotropy fluctuations in the range of a few per- 
cent are not unreasonable in nm-scale devices. Calcu- 
lations for transition-metal clusters show that single spin 
flips can produce a significant change in the magnetic 
anisotropy energy jl5| , and measurements of Gd clusters 
indicate that anisotropy energies can vary significantly 
for clusters differing only by a sing le atom Q. 

We now turn to the low-energy excitations observed 
in (PI). It is natural to ask |],[| whether these 

might correspond to spin wave transitions of the form 
\si,Si) N — » \s f,mf) N±1 for different rrif values. How- 
ever this does not seem to be the case, for three rea- 
sons: (i) It can be shown that only two transitions 
(namely S07S0} — > \sq ± 1/2, so ± 1/2)) have significant 
weight ■ Resonances associated with final states 
|s/,toj) that differ only in m/ would (ii) have a spacing 
of order [~ 0.01 me V], which is significantly smaller 
than observed, and would (iii) exhibit a systematic in- 
crease in the magnitude of their slope (oc \si — mf\) for 
high magnetic fields that was not observed in experiment. 

Nonequilibrium. — Since the large density of resonances 
(PI) cannot be explained by equilibrium transitions (nei- 
ther single particle excitations nor spin wave excitations) , 
we must explore nonequilibrium effects: In general, N- 
electron states other than the ground state can be pop- 
ulated during the process of current flow, and this may 
affect the experimental tunneling spectrum [p^p|. Fig- 
ure H illustrates the consequences as applied to a ferro- 
magnetic grain. Even if a first tunneling event causes a 
"charging" transition from the A^-electron ground state 
\G) N to the (A r ±l)-electron ground state |G) Ar±1 , it may 
be energetically possible for the subsequent "discharging" 
tunneling transition to return the particle to an excited 
TV-electron state \a) N instead of IG)^, provided the ap- 
plied voltage is sufficiently large, eV > E^ — Eq . Like- 
wise, further charging and discharging transitions may 
allow any of a large ensemble of states to be occupied at 
higher and higher levels of an energy ladder, terminating 
only when an energy-increasing transition requires more 
energy than the applied voltage provides. As the volt- 
age is increased, the total current (conductance) may in- 



crease stepwise (show peaks) when thresholds are crossed 
to allow higher-energy transitions up the nonequilibrium 
ladder, thereby changing the occupation probabilities of 
the ensemble of nonequilibrium states and opening new 
tunneling channels. 

In a ferromagnetic particle, in addition to the nonequi- 
librium occupation of single-electron states discussed pre- 
viously for nonmagnetic particles JlS] , nonequilibrium 
spin excitations are possible, too, if the spin-flip rate T s f 
is smaller than the tunneling rate r tU n jl9| ■ Fig.||A illus- 
trates a ladder of transitions between ground states \s,s) 
p0| of multiplcts with different total spin s, leading to 
spin accumulation on the grain. Fig |^B shows the cor- 
responding differential conductance, calculated by solv- 
ing a master equation for the population of the states 
of the ladder. The resonance peak spacing (SE ICS ) and 
the number of peaks (n res ) for such a ladder can read- 
ily be calculated (using Eq.(2) of 0), and depend on 
whether the charging transition adds/removes an elec- 
tron to/from the grain (to be distinguished by an index 
p = ±1), and on whether it is a majority/minority elec- 
tron (to be distinguished by an index a = ±1). One 
finds SE rcs — d mm — U/2 for pa — —1, and rf ma j — U/2 
for pa = 1. Using the model parameters estimated 
above with sq ~ 1000, the first quantity gives a spac- 
ing of ~ 0.2 meV, as is observed. The second quantity 
is larger, r* 3.6 meV A detailed analysis jlTj shows 
that n ros — 1 equals the smallest integer larger or equal to 

2i% hrcsh ~ 2B„Aiffi ot 

(Bp — B- P )(d ma j + d min — U) — a(d ma j — d m i n ) 

where B p = [1 + C p /C_ p ] _1 contains the ratio of junc- 
tion capacitances involved in processes p and —p, AE^ ot 
is the energy difference between states \sq + ap, sq + ap) 
and |so,So)j and E^ rcsh is the threshold charging en- 
ergy (energy of the first peak in the differential conduc- 
tance). The prediction that rt rcs increases linearly with 
£^ hresh @ is in qualitative agreement with Fig. 2 of [|, 
and could be checked explicitly in future devices with 
gate electrodes, which would allow £^ rosh to be tuned. 

A nonequilibrium scenario can also account naturally 
for the fact (P3) that the vast majority of the observed 
transitions within a given sample shift in energy with a 
similar slope for large magnetic fields. This will happen 
when all the nonequilibrium threshold transitions corre- 
spond to tunneling events with the same change of S z 
(see Fig. ||), and therefore the same Zeeman shift. 

In conclusion, we have proposed a phenomcnological 
model for nanoscale magnets that treats magnetic inter- 
actions within a many-electron picture. Its parameters 
were estimated from bulk properties of Co or experiment, 
except for the total spin s and the strength of anisotropy 
fluctuations, which were used as free parameters. The 
model offers a framework for understanding recent ex- 
periments measuring the discrete excitations of magnetic 
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nanograins, provided that we assume (i) anisotropy fluc- 
tuations of a few percent between different eigenstates 
within the same nanograin, and (ii) nonequilibrium spin 
accumulation. 
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FIG. 1. Tunneling energies A-Ei for 7iz co + Hum, plotted 
as functions of /i/fcjv sweeping positive to negative, illustrat- 
ing the effects of anisotropy fluctuations (Sk± = fcjv±i — fcjv) 
for the transitions from Sj) to (A) |st 
and (B) k + §,Sj + %) N+1 . 
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FIG. 2. Nonequilibrium spin accumulation in a ferromag- 
netic nanoparticle: (A) Tunneling transitions can cause an en- 
ergy ladder of states with different total spins to be populated 
(so denotes the ground state of the spin-so multiplet, etc.). 
(B) The corresponding differential conductance as function 
of energy, normalized by its first maximum and calculated by 
standard methods [^. (Parameter choices: h = 0, 7i a ms = 0, 
so = 1000, AE^ tot = 0, T = 80mK, B p = 0.4, a tunnel junc- 
tion resistance ratio of Rl / Rr = 0.1, and a Coulomb blockade 
threshold of 7 meV (see sample 3 in Fig. 2 in we neglected 
energy and spin relaxation, and assumed that charging tran- 
sitions add minority electrons, so that SE las = d m m ~ U/2.) 
No significance should be attached to peak heights here, since 
they depend on (unknown) tunneling matrix elements, which 
for simplicity we took all equal. 



